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Dynamics of decoherence: universal scaling of the decoherence factor 
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We study the time dependence of the decoherence factor (DF) of a qubit globally coupled to an 
environmental spin system (ESS) which is driven across the quantum critical point (QCP) by varying 
a parameter of its Hamiltonian in time f as 1 — t/r or —tfr, to which the qubit is coupled starting 
at the time t —>■ —oo; here, r denotes the inverse quenching rate. In the limit of weak coupling, we 
analyze the time evolution of the DF in the vicinity of the QCP (chosen to be at t = 0) and define 
three quantities, namely, the generalized fidelity susceptibility (defined right at the QCP), 

and the decay constants ai(r) and 02 ( 1 ") which dictate the decay of the DF at a small but finite 
t(> 0). Using a dimensional analysis argument based on the Kibble-Zurek healing length, we show 
that Xf{f) as well as Qi(r) and 02 (t) indeed satisfy universal power-law scaling relations with r 
and the exponents are solely determined by the spatial dimensionality of the ESS and the exponents 
associated with its QCP. Remarkably, using the numerical t-DMRG method, these scaling relations 
are shown to be valid in both the situations when the ESS is integrable and non-integrable and also 
for both linear and non-linear variation of the parameter. Furthermore, when an integrable ESS is 
quenched far away from the QCP, there is a predominant Gaussian decay of the DF with a decay 
constant which also satisfies a universal scaling relation. 

PACS numbers: 


In the context of quantum computation and informa¬ 
tion mm, one of the major issues is the study of deco¬ 
herence mm, namely, the loss of coherence in a quan¬ 
tum system due to its interaction with the environment. 
To investigate the environment induced decoherence of a 
qubit in the vicinity of a quantum critical point (QCP) 
mm of the environment, a paradigmatic model known as 
the central spin model (CSM) [ 3 | 8 ] has been generalized 
to the context of a quantum phase transition [^. 

In the CSM a central spin (CS) or a qubit is glob¬ 
ally coupled to an environmental quantum many body 
system, usually chosen to be a quantum spin system, re¬ 
ferred to as the environmental spin system (ESS) in the 
subsequent discussions. The ESS is initially in its ground 
state while the CS is in a pure state; the global coupling 
between the qubit and the environment is so chosen that 
the subsequent time evolution of the initial ground state 
wave function of the ESS occurs along two channels dic¬ 
tated by two different Hamiltonians. Even though the 
qubit is initially in a pure state, it has been shown that 
it loses its purity (almost completely) [5j when the ESS 
is close to its quantum critical point (QCP). Question we 
address in this letter is as follows: what happens when 
the ESS is slowly driven across its QCP? Is there a uni¬ 
versality associated with the dynamically generated de¬ 
coherence quantified by the decoherence factor (DP) of 
the CS especially in a limit when the coupling between 
the CS and the spin chain is weak? 

Although we shall consider more generic non- 
integrable models in this letter, let us first illustrate the 
basic idea using the transverse XY chain consisting of N 
spins m as the environmental spin chain; the model is 


described by the Hamiltonian 

N 

hW = -X (1) 

where cr’s are the standard Pauli matrices, Jx{Jy) is 
the ferromagnetic nearest neighbor interactions along the 
x{y) directions and h is the transverse field; this spin 
chain is coupled to the spin- 1/2 qubit by a Hamiltonian 
Hse- In the following, we shall set + Jy = 1 and the 
anisotropy in interaction i.e., Jx — Jy will be denoted by 
the parameter 7. The phase diagram and different phase 
transitions of the model Q are presented in the Pig. Q. 

Let us first assume that the qubit is coupled to the tun¬ 
able transverse held of Eq. 0 through the Hamiltonian 
Hse = ^i^s, where af is the i—th spin of the 

XY chain and ag represents that of the qubit with 5 being 
the coupling strength. [Notably, whenever the qubit is in 
the up (down) state, the held h of the ESS gets altered 
io h — 5 (/i -I- 5 ).] We choose the qubit to be initially (at 
t —>■ —00) in a pure state \(j)s{t —t —00)) = ci| f) + C2I i) 
with |ci/ -I- IC2/ = 1 , where | f) and | |) represent up 
and down states of the CS, respectively, and the environ¬ 
ment is in the ground state \ 4 >Eit —t —00)) = \ 4 >g)- The 
initial state of the composite Hamiltonian He + Hse, 
at t —>■ —00, is then given by the direct product \tj;{t 
—00)) = \ 4 >s{t —t —00)) C) |/)g). It can be shown that at 
a later time t, the composite wave function is given by 
\ip{t)) = Ci| t)C)|</>-s(Q) + C2| i) 0 |/_(t)), where |/±) are 
the wavefunctions evolving with the environment Hamil¬ 
tonian H^'^{hzLS) given by the Schrodinger equation 
id/dt\(j)±) = H^{h ± ( 5 )|</>±). We therefore hnd that the 
coupling 5 essentially provides two channels of evolution 
of the environmental wave function dictated by two trans- 
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FIG. 1: (Color online) The schematic phase diagram of the 
transverse XY chain given in Eq. Q in the {h — 7) plane with 
different quenching paths denoted by arrows. The Ising tran¬ 
sition line at /i = ±1 for arbitrary 7 and the anisotropic tran¬ 
sition line at 7 = 0 (with \h\ < 1) meet at two mulcti-critical 
points (MCPs) h — ±1 and 7 = 0 . When the transverse field 
h = 1 — tjr is quenched (e.g., setting 7 = 1 ) with t from a 
large positive value to close to the QCP (t > 0), one finds 
Xf{t) ~ , ai(r) ~ r°, and 0 : 2 ( 1 ") ~ in the limit 

of small (5; identical scaling relations are obtained when 7 is 
quenched across the anisotropic critical line. Finally when 7 
is quenched linearly across the MCP with h = 1, one finds 
Xf{t) ~ r^/^,oi(r) ~ and 02 ~ As shown in the 

text and in the supplementary material that all these scal¬ 
ing relations are predicted by universal scaling form given in 
Eq. 

verse XY Hamiltonians with the transverse field h+6 and 
h — S, respectively. It is straightforward to show that the 
DF defined through the relation D{t) = \{(j)+{t)\(j)-{t))\‘^ 
measures the purity of the reduced density matrix of the 
qubit; any non-zero value of In D(t) implies that the qubit 
is in a mixed state. When the transverse field is close to 
the quantum critical value there is a sharp dip in D{t) 
which on the one hand, detects the existence of a QCP of 
the ESS and on the other, establishes that the qubit loses 
its initial purity almost completely in its vicinity |^. 

The quantity D{t) is also known as the Loschmidt echo 
which has been studied in recent years (both at zero [5] 
and finite temperatures El) in the context of decoher¬ 
ence in equilibrium [SIIMI] and non-equilibrium situ¬ 
ations and is closely connected to the dynamical 

phase transition nnun], the statistics of work done [12] 
and the entropy generation in a quench P51 - E5] . 

We assume a weak coupling between the qubit and 
the environmental Hamiltonian; furthermore the qubit 
is always globally coupled to time-dependent part of the 
driven Hamiltonian. Under these circumstances, we now 
assume the generic situation when a parameter of the 
ESS Hamiltonian is quenched A = —tjr, with t starting 
from a large negative value; here, A is the deviation from 



FIG. 2: Time evolution of the decoherence factor for the 
transverse Ising chain (TIG) as obtained using the t-DMRG 
method. The transverse field is driven as h = 1 — t/r. The 
scaling of Xf{t), 01 ( 1 ") and 02 ( 1 ") match nicely with the the¬ 
oretically predicted scaling relations given in Eq. §. 

the QCP (e.g., when the field h in model 0 is quenched, 

A = h — 1), and r(^ 1) is the inverse quenching rate. 
Consequently, one has two channels of evolutions of the 
initial ground state of the ESS dictated by two Hamilto¬ 
nians with parameters A -I- d and X — S, respectively. 

In general, this quantity D{t) can be expanded for suf¬ 
ficiently small S and t as 

D{t) « 1 - (^xf{t) + ai{T)t + (2) 

or equivalently, 

^ lnD{t) « - I^XP’(r) -b ai(r)t -b ia 2 (T)t^^ <5^ (3) 
with 

1 Id™ 

XFir) = -^lnD{0), a:„(r) = -^—(InD(t))|t^o, 

where L is the linear size of the system. Eigurej^ shows 
the time evolution of D{t) in the transverse Ising chain 
(TIC) derived from Eq. Q with 7 = 1 (i.e., Jy = 0) 
driven as h = 1 — t/r. One can see that D{t) follows 
Eq. (§ with xf and am changing with r. It should 
be noted that there is a prominent Gaussian early-time 
decay of D{t) (i.e., a\ = 0) as predicted in the Ref. 

[26] only when II(t) is measured between the ground 
state of unperturbed Hamiltonian; in the present case 
|^+(t = 0)) ^ = 0)), and hence we indeed have 

a linear term in Eq. ^ which is also numerically estab¬ 
lished below and analytically established in the supple¬ 
mentary material (SM). 

In this paper, we present the universal scaling of 
these quantities namely, the quantities Xf(t), Q:i('r) 
and 02 (t). We note that lnT>(A = t = 0) = 
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ln(|((^_i_(t = 0)|(/)_(i = 0))p) in fact defines the fidelity 
m between two states evolved from the same initial 
state along two different channels up to the value A = 0. 
In the limit i5 —>■ 0, lnI?(A = t = 0) can be expanded in 
the form lnZ3(A = t = 0) = where Xf{t) 

can be viewed as a generalized fidelity susceptibility (see 
the supplementary material (SM)). 

In order to resolve the question raised above concern¬ 
ing the scaling of Xf{f), and 02 (t), we resort to 

dimensional analysis. When the ESS is slowly ramped 
across its QCP, according to the Kibble-Zurek Scaling 
(KZS) argument, the non-adiabatic effects be¬ 

come prominent only in the vicinity of the QCP; this is 
a consequence of the diverging relaxation (healing) time 
at the QCP. As a result, there is a characteristic time 
scale of the problem i ~ ^ jg ^j^g corre¬ 

lation length exponent and 2 is the dynamical exponent 
associated with the QCP across which the ESS is driven. 
This follows from the simple adiabatic impulse argument 
which relies on the assumption that at t = t, the rate of 
change of the Hamiltonian equals the characteristic (heal¬ 
ing) time of the problem. So far as low energy modes are 
concerned, one can assume that within the time window 
—t<t<t, known as the impulse region, the wave func¬ 
tion remains frozen to the adiabatically evolved state at 
t = —t, resulting in defects in the final state. One there¬ 
fore finds a characteristic length scale, known as the heal¬ 
ing length L that scales as L ~ assuming there 

is one defect within the healing length L, one immediately 
finds the KZS of the defect density n in the final state 
following a slow ramp given by n ~ lfL‘^ ^-vd/{vz+\) ^ 
We note that the proposed KZS form has been verified 
and modified in several situations [5T1 - HD] (For review see 
[inumiia). We shall demonstrate below how the scaling 
of L leads us to the universal scaling of xf{f), ch(t) and 
02 (t) that we want to establish. 

To derive the scaling, we first use the definition of the 
critical exponents v and z to find the scaling dimensions 
t ^ and 6 (which is equivalent to A) ^ . De¬ 

manding that lnT)(A = t = 0) must be dimensionless, we 
find 


Xf{t) ~ 

ar„('r) ~ (m=l,2). (4) 

Remarkably, the exponents appearing in scaling relations 
are entirely determined by the universal critical expo¬ 
nents and the spatial dimension d. As shown below, these 
relations hold true for all integrable as well as the non- 
integrable situations discussed here. Remarkably, the 
scaling of Xf{t) also follows from the well established 
scaling of the fidelity susceptibility associated with the 
ground state fidelity (see the SM). 

Question may arise what would happen if one consid¬ 
ers a non-linear drive X{t) = — |t/T|''sign(t) of the ESS 


where sign(t) denotes the sign function. In this case, L 
scales as [37j[38] and following a similar line 

of arguments, one finds 


XFir) ~ ^R2-d.)/(r..+ l)^ 

amir) ~ ^r{2-du-mz^}/irz^+l} (^= 1 ^ 2 ). ( 5 ) 

Equations Q and ([^ contain the central result of our 
paper. 



(a) 


(b) 


FIG. 3: (Color online) (a) Scaling of xf, cei, and a .2 with re¬ 
spect to T for the TIC (Eq. 0 with Jy = 0) with a non-linear 
ramp (r = 2) of ft. We find an very good agreement with an¬ 
alytically predicted scaling in Eq. Q. (b) Scalings of Xi, Oi 
and 02 for the quadratic ramp (r = 2) of the longitudinal 
field in the non-integrable TIC given in ([^ with ft = 1; the 
numerical results as obtained using t-DMRG are in excellent 
agreement with the predictions as in Eq. (101. 


Let us investigate above scaling relations for the inte¬ 
grable model 0 when the transverse field ft is quenched 
as ft(<) = 1 — t/r with Jy = Q and the spin chain is 
quenched across (but close to) the QCP at ft = 1 (i.e., at 
t = 0) with exponent z = f = d = 1: Eqs. (0 and 0 
predict: 

XF(T)~ri/^ ai(T)-T°, a2('r)~T"^/^ (6) 



t 


FIG. 4: Early-time evolution of the decoherence factor for 
the non-integrable TIC in the longitudinal field 0. The 
transverse field is fixed at ft = 1 and the longitudinal field 
is ramped as fti = —t/r. Inset shows scaling of xl, cti and 
02 - Results are in good agreement with the prediction 0. 
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for the linear ramp, while 

Xf(t) ~ ai(T)-T°, a2(r) ~ (7) 

for the non-linear ramp. The scaling relations in Eq. ®> 
has been verified numerically (see Fig. (H) and estab¬ 
lished using exact analytical calculations as shown in the 
SM. Scaling relations for the non-linear ramp presented 
in Equation Q is also numerically verified with an ex¬ 
cellent agreement as shown in Fig. [^a). 

We now proceed to examine the universal scaling when 
ESS is non-integrable. We consider the situation when 
we have a ferromagnetic Ising chain in both transverse 
and longitudinal fields [i5] : 


Hh = -Yl E - E (8) 

i i i 

clearly the Hamiltonian is at the quantum critical 
point when = 0 when it represents a critical quantum 
Ising Hamiltonian obtained by setting h = = 1 and 

Jy = 0 in Eq. Q. The dynamical exponent z = 1 and 
a gap opens up in the spectrum for small non-zero 
as , leading to a value of f = 8/15. Considering 
a similar coupling with the qubit of the form HgE = 
—5 with i5 —)■ 0, we study the quenching of the 

ESS close the integrable quantum critical point (/il = 0) 
by varying = —^It starting from a large positive value 
of Hl- Scaling relations as obtained from Eqs. Q and 
@ are 

xAt) ~ ^ ~ (9) 


for the linear ramp and 


Xf{t) 
0^2 (t) 


22r/(8r-|-15) 

' ; 

6r/(8r-|-15) 


ai(r) 


^14r/(8r+15) 


( 10 ) 


for the non-linear ramp /iL(t) = —|t/T|’'sign(t). Figure [4| 
shows the early-time behavior of the DF for the model (1^ 
following a linear quenching of as obtained using the t- 
DMRG method. The scalings of xf, cti, and q ;2 extracted 
close to t = 0 are good agreement with the prediction in 
Eq. (§. Figure[§b), on the other hand, shows the scaling 
of xf, Oil, and q ;2 for the non-linear ramp; numerical 


results and theoretical prediction in Eq. (101 with r = 2 
are in excellent agreement. 

Let us now address the issue that what happens to 
the DF when the system is quenched far away from the 
QCP, i.e., h in model 0 is changed from a large nega¬ 
tive to a large positive value. For an integrable model re¬ 
ducible to decoupled two-level problems (like the model 
Q or the Kitaev honeycomb model [44]), one can an¬ 
alytically establish that in the limit f 2> 0, InZ? ^ 
i.e., there is a prominent Gaussian de¬ 
cay with time niiiT]. Using the dimensional analysis ar¬ 
gument presented above, we immediately conclude that 


52(t) ~ ^{ 2 -dv- 2 zv)/(zi^+i) linear quenching of A 

while for a non-linear ramp 52 ('r) ~ j-r{ 2 -diy- 2 zi/)/{rziy+i)_ 
It is noteworthy that whenever vz = 1, the scaling of d 2 
is the same as that of n for both r = 1 and r ^ 1; oth¬ 
erwise, the scaling of n and 52 are completely different 
m The universal scaling formula of 52 (f) explains all 
the integrable situations discussed in earlier studies. (See 
the SM for relevant discussion.) 

In summary, using the GSM and working in the weak 
coupling limit, we have provided a universal scaling re¬ 
lation of the generalized fidelity susceptibility Xf{t) as 
well as the decay constants Q!i(t) and a 2 (T) with the 
inverse rate r in the vicinity of the QGP across which 
the ESS is driven. We have used a dimensional analy¬ 
sis arguments based on the healing length L. Remark¬ 
ably, the exponents are solely determined by the spatial 
dimension of the ESS and the associated critical expo¬ 
nents. These scaling relations are verified for both in¬ 
tegrable and non-integrable models as well as for linear 
and non-linear ramps. Additionally, we have also derived 
a universal scaling relation for the decay constant 52 (t) 
describing the Gaussian decay of the DF with time when 
an integrable ESS is quenched far away from the QGP. 
We reiterate that the DF is equivalent to the Loschmidt 
echo and the latter has been experimentally observed in 
the TIG using NMR technique [ISj. Therefore an ex¬ 
perimental test of our scaling theory on the DF can be 
realized in a large scale NMR quantum simulator. Fur¬ 
thermore, Kibble-Zurek scaling is already being explored 
in optical lattices [15]. 
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nancial support. 
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In this supplementary material, we shall illustrate how to justify the expansion of In D{t) close to t = 0 and derive the scaling 
of Xf (''")) Qi(''') 8-nd a 2 (r) for a transverse Ising chain using scaling arguments and exact analytical calculations. 

In order to estimate the DF in the limit t —>■ O-f, we first define the generalized time-dependent fidelity susceptibility XF{t) 
as 


= (SI) 

and expend this quantity into Taylor series of t as follows: 

XF^t) = Xf{t) -j ai(r)t -|- ia2(r)t^ -!-•••. 


(S2) 
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Therefore, the DF in the limit t —>■ 0+ comes out to be 

D{t) = ~l-3^LxF{t) 

1 


Xf{t) + ai{T)t + -a2{r)t 


In the following, we shall derive the scaling of Xf{f), cii('r) and a 2 ('r) in integrable situations. 


(S3) 


THE SCALING OF ai(r) AND a 2 {T) FOR INTEGRABLE SITUATIONS 


We shall first derive the scaling of XF('r) or the generalized fidelity susceptibility using the established scaling of the 
ground state quantum fidelity. As discussed in the main text, let us assume a linear quenching scheme A = —t/r (for 
the quenching of the transverse field h in the transverse Ising model, X = h — 1) which stops at the critical point at 
X = t = 0. We then calculate the logarithm of the overlap or fidelity, lniA(A = t = 0) = ln(|((^_|_(t = 0)|(()_(t = 0))|); 
expanding up to the second order in S, one finds ln(|(0+(t = 0)\cj)_{t = 0))|) ^ —S‘^L'^xf{t)- It is well established 
that at the QCP the ground state quantum fidelity susceptibility scales as where v is the correlation 

length exponent, d is the spatial dimension of the system, and L is the linear size of the system. In the present 
problem of slow ramp, the characteristic length scale, as discussed in the main text, is given by the healing length 
L Using L in the scaling of the ground state fidelity susceptibility one arrives at the scaling relation 

Xf{f) j{'^-vd)l{vzFi)^ which has been obtained in the main text using a dimensional analysis approach. We note 
that for this scaling to be valid, one must have ^ L ^ L. 

We shall illustrate the above scaling considering the environmental Hamiltonian to be a transverse XY chain 


He 


E 


f 7 


i+l 


+ Jycrfcrf+i + h{t)a^) , 


(S4) 


and we first consider the Ising limit obtained by setting Jy = 0; we also choose Jx = I- The coupling Hamiltonian 
between the system and the environment is given by 


L 

Hse = -(5^ erf erg, 
2=1 


(S5) 


where erf is the i-spin of the environment chain and ag represents the spin of the system. Hereafter we assume a 
linear ramp of the transverse field: 

X{t) = h(t) — 1 = -, (S6) 

T 

where r stands for the inverse rate of a ramp. 

According to the exact solution presented in Ref. [32] for L —>■ oo, the decoherence factor is given by 

D{t) « exp J dk\nF-^{t)^ , (S7) 


Fk{t) = \u+*{t)u^ (t) + v+*{t)v^ (t)p, (S8) 

where and are the solutions of the time-dependent Bogoliubov-de Gennes equations 

= 2{h{t) ± 6 — cos k)u'^ + 2 sin k , (S9) 

= 2 sin k — 2{h{t) ± 6 — cos k)v'^. (SIO) 


Using the notation, t(j_ = —4Tsink{h{t) ± J — cosk), t' = drsin^fc, z± = t'^e^'^Hj ^nd n = —ir' 14, the exact 
solution of the above equations under the initial conditions it^(—oo) = 1 and t>^(—oo) = 0 is obtained as 

W±(t) = e-’^^'/iegrrrU 

X [(1 -I- n)'D^ri- 2 {iz±) F iz±'D^n-i{iz ±)\, 


(Sll) 
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(S12) 


where Vmiiz) is the Weber’s parabolic cylinder function. We note that Fk{t) = 1 when (5 = 0, because Fk{t) reduces 
to the norm of the wavefunction. 

Hereafter we fix t = 0 so that h{t) = 1 and assume r ^ 1. Under this assumption, one can show that only the long 
wave lengths with k ^ <C 1 give contribution to the integral of InU^T^ in Eq. (S7). Therefore, in the present 

4t/c2, 


situation, one has « — 4t/c(±^ + k^l2), t' 


z± « —2-y/T(±5 — fc^/2), and n « —irk 


Let us consider the situation where b ^ ^ 1. Noting that and depend only on r', z± and n, one can 

see that InFjr^(O) is given as a function of r^^'^k, and Therefore the Taylor expansion of InF^^(O) up 

to the order of tS'^ leads to InFjr^(O) « (^(r^^^fc, where (j){x,y) is a function of x and y. We remark here 

that the linear term of (5 vanishes since Fk is an even function of S. To carry out the integral of InE^^, one can safely 
shifts the upper limit to infinity because there is no contribution to the integral from k ^ . Consequently, one 

obtains 


--ln£i(A = t = 0 ) - / dfc(/)(r^/^fc,T^/ 2 /c^)r (52 

^ Jo 

POO 

^ .j.i/2^2 / dk4)(k,T~^^‘^k^) 

Jo 

^ .,- 1 / 2(52 |y dk (j){k,0) + O 

- (S13) 

which gives Xf(t) = — lnZ?(A = t = Q)/[L 6 '^) ^ Noting that the critical exponent v associated with the QCP 

at /i = 1 is 12 = 1 and spatial dimensionality d = 1, this scaling is in perfect congruence with the expected universal 
scaling behavior of xf- In non-linear ramping case, A = — |t/T|’'sign(t), we can proceed using a similar arguments 
and establish that Xf{t) ^ ^ 

Now, we shall calculate Fk{t) explicitly by expanding the Weber function, Vniz), near z = 0. One can show that 
2 <C 1 demands <C VF. We work in the limit of 6‘^t <C 1. The power series expansion of 'Dn{z) for small 2 is given 

by 


Fniz) 



2"2^y^z 2 ^~^2^7r(l + 2n)z^ 

r(^) ■ r(i-§) 


The solution of time-dependent Bogoliubov-de Gennes equations for t = e are given by 




g-7i-r7l6gi7i-/4 


■ 2 S07 

.r(i + f) 


t'z±^Jtt 2 2 ^ (1 — 2 n ) z'^^/tt 2 2 

4^2 r(i + §) ^ ^ r(i + t) 


v^{t) ~ 


/I6 9 - 


1 

•\/2 iz± 

(1 - 1 - 2 n) z 4 

2 ^ 

V 2 

[r(i + f) 

r(i + f) 

4 r(l + f)_ 


(S14) 


(S15) 


(S16) 


where z± is defined by Using the above expressions one can find out the quantity u'l*{t)uf. (t) + 

which is given by 






Z^ ai + Z- 02 + Z^Z- 03 + (z^)^ 04 + z'^ 05+ z+zi 06 -I- (z!j_)^z_ 07 


+ 0(4), 


(S17) 
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where 


8r(i + ^f)r(i + f) 8r(i + f)r(i + ^f)’ 
02 =0^, 

ttt' ttt'^ 

^ 4 |r(i + f )|2 + 32 |r(i + f)| 2 ’ 

7rT'(l + 2 n*) 7r(l — 2 n*) 

="32 |r(i + f )|2 + |r(i + f)| 2 ’ 


0-5 — ^4 5 


iTTT'{l + 2n)-\/2 


7 rT'(l — 2 n) 


32 r(i + )r(i + f) 4 v^ r(i + f )r(i + f) ’ 


07 = Og . 


Therefore, Tfc(f) is given by Tfc(f) = 1 + + C'fe(f) + -Dfc(f), where 


7lfc(f) = 2e“’"^'/® |Re(z+ai)| + |Re(z_a2)| , 

Bk{t) = 2e“’"'^|Re(z+z_a3)| + |Re(z+04)1 + |Re(z^05)1 + 2|Re(z_z+a20*)| + |z+p|aip + |z_p|a2p , 

Cfc(f) = 2e“'^'^|Re(zi^z^a6)| + |Re(z^z* Oy)| + 2e“'^'^ |Re(|z+p2_a4a3)| + |Re(|z_pz+a 203 )| 


+ |Re(z+z+aia4*)| + |Re(z+(z*) aiOg)! + |Re(2_z+a204)| + |Re(z_(z_*) 0205)^ , 
Dk{t) = e-""'/4[|^^|4|„^|2^|^_|4|^^|2^|^^|2|^_|2|^^|2l + 2e-^^'/^\\Re{izlf (zlf a,4)\ 


+ |Re(|z+p(z* )^aia6)| + \Re{z+z*_\z+\^ai4)\ + |Re(z+z* |z_pa20g)| + |Re(z|zl|2_^0305)! 

+ |Re(z+| 2 _p 0207 )| + |Re( 2 +z_| 2 :+pa 304 )| . 

The necessary algebra leads us to the conclusion that Ak{t) = 0, which establishes the fact that there is no linear 6 
term present in Fk(t). One can get the foilwing contribution from Bk(t) 


Bk{t) ~ -S^ [2TT'e-’^^'/2 + 

L J L ^ ^ 


Similarly, Ck{t) gives 0{z^) term, 
Ck{t) ~ -6^ t 4l4-^-r'/s 


3 ^-3/2 -^t'18 




967 rT '2 967 rT' 1 40 S^nt 

64 |r(i + f)r(i + f)| ■ 8 |r(i + f)r(i + f)|J ” 3 |r(i + f )r(i + f)| 

967 rr '2 Ottt' 1 


64|r(i + ^^)r(i + f)| |r(i + f)r(i + i^)|J |r(i + f)r(i + ^ 


The fourth order term in z± is given by Dk{t), 


Dk{t) ~-5‘^t^ e ^Qfi ++ e x ^a2 +/?2r'^^ - — [^3 + Vi'^'+ 

+ — ^^4 + 772/ + /34'r'^^ - ^as + ri^r' + ^ag + ?? 4 'r' + / 3 gr'^^V) 


as + Vit' + /Jst'^ 




where a’s, /3’s and ry’s are the numbers coming from the T-function. 
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Using the approximation of ln(l + cc) ~ a; for a; <C 1 and expanding near the critical point h = 1 and the critical 
mode fc = 0, we can now compute the DF by evaluating the following integral over momentum space: 


/“ dklnFkit ) ~ - 5 '^ + jdt + _ <^^2^-372 _ ^^3^-2 _ ^^4^-572 _ 5,^4^372 

~ -(5^ (q!T^/^ + /3t + H- 


(S22) 


In the last line we have included only the leading order contributions to the integral. Here, a, j3, ■ ■ ■ are the numerical 
prefactors coming from the integartion of the dimensionless F-functions and the dimensionless exponentials. 

The DF in the limit t = e, e —>■ +0, is therefore given by 


D{t) 


exp 


~ exp - 


-I 

27r Jo 

6^L 


Ito 


~ 1 - 


27r 


dk lnFfc(t)J 
+ /3t + 
+ /3t + 


.172 


ai(T) 


and 02 (r) 


.-172 


(S23) 


as predicted by dimensional analysis argument 


Clearly, one finds that x_f(t) 
given in Eq. (§. 

Referring to the Fig. 1 of the main text, we consider now the case when 7 = J^, — Jy is ramped as 7 = —t/r with 
Jx + Jy = h = 1 fixed in Eq. (S4). In this case, one has a ramp along the gapless Ising critical line across the MCP. 
Using the notion of the dominant critical point and exponents d = v = 1 and z = 2[gi|l], Eq. 0 yields 


XF(r)--r/'^, ai(T)-T 02(4-) 

for the linear ramp (see Fig. [sT]^a) ) while 


^-1 


Xf(7 


r7(2r+l) 


ai(T) 


^-r/{2r+l) 


02 ( 4 ") 


^-3r7(2r+l) 


(S24) 


(S25) 


for the non-linear ramp (see Fig. |Sl[ b)) . 

We shall verify the above scaling form of Xf(t), and a 2 (r) by explicitly calculating the DF. The interaction 

Hamiltonian between the environment and the qubit is accordingly modified as Hse = —{5/2) Y/,i{^i^i+i—'^l^l+i)^s- 
The coupling <5 therefore provides two channels of the temporal evolution of the environmental ground state with 
anisotropy 7 ^ and 7 — d. The time-dependent Bogoliubov-de Gennes equation governing the time evolution of 

environmental state through the gapless line are given by 


= 2 ( 7 ( 7 ) =F (5) sin/c wj-I-2(/i - cosfc) 
= 2(/i — cos fc) — 2 ( 7 ( 7 ) =F (5) sinfc 


(526) 

(527) 


Accordingly, we use the notation, 7(j_ = 4t(7i — cos k){^{t) =F 5), t' = 4r(h — cosfc)^/sinfc, z± = ^'^/y/F and 

n = —ir'/A. 

The solutions are the same as obtained for the case of transverse field quenching given in Eq. (Sll) and Eq. (S121. 
In order to estimate Fk{t) with 7 = e, e —>■ -1-0, one has to use the following approximation: 7(j_ <C Vt' and S'^rk <C 1. 
Following the same line of algebra one can find out a similar expression of ul*{t)u'^{t) + {t)v//[t) with a modified 

^gapless _ 2-v/rsinfc(7/r =F <5) exp(7F/4). One can note that for the Ising model zj™® = 2y^(7/T =F <5) exp(7F/4). 
As a result, = Vsinfc z^™®. In this case, we find Au{t) = 0. The O(z^) contribution is given by = 

sinfc "®. Similarly, Cf = sin^/^ fc cf"® and Df?'®- = sin" fc 

One can therefore calculate the integral of Fk (7) by using the logarithmic approximation and expanding momentum 
near the critical mode, sin k cs k and t' ~ fc®r 


c7fclnFfe(7) ~ -5" {ar^/^ + F 


r 

JO 

~ —(5" -!-•••. 

The DF D(7) in the vicinity of MCP at 7 = e, e —>■ -1-0 is given by 

n(t) ~ 1 - ^ (ari/3 + /37r-i/® + 7y7"r-i). 

ZTT 


(S28) 


(S29) 
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FIG. SI: (a) Early-time evolution of the decoherence factor for quenching the parameter 7 as —I/t of the Hamiltonian 0 
with the transverse field is fixed at h = 1; the left panel shows scaling of xl, oti and a 2 - Results are in good agreement with 
the prediction in Eq. (S24l. (b) Numerical results confirm the scaling for the non-linear case given in Eq. (S25l. 


This expression of D{t) clearly verifies the scaling of Xf(t), and 0 : 2 (t) predicted above in Eq. (S24). 

On the other hand, for an anisotropic quenching when 7 is quenched with \h\ < 1 and the gap vanishes for the 
critical mode = cos“^ /i, we have t' ~ 4t\/1 — and sinfc ~ VT^-T?; using these parametrizations, one can 
retrieve scaling Xf{t), <ai('r) and 0:2 ('t) similar to those obtained in the /i-quenching case. 


INTEGRABLE MODELS: SCALING OF In D{t) FAR-AWAY FROM THE CRITICAL POINT {t > 0) 

Let us now address the issue that what happens to the DF when the system is quenched far away from the QCP, 
i.e., h in model 0 is changed from a large negative to a large positive value. For an integrable model reducible to 
decoupled two-level problems (like the model 0 or the Kitaev honeycomb model [3S]), one can analytically establish 
that in the limit t 0, Inl? ~ 5^L'^a 2 {T)), i.e., there is a prominent Gaussian decay with time [321 [ 6 ]. Damski 

et al. proposed that for an integrable ESS, Q; 2 ('r) is likely to scale as n as determined by the KZS relation. In 
a subsequent work. Nag et al. [3S] showed that this conjecture is not necessarily true even for the above integrable 
model when the system is quenched along a gapless critical line of the model 0 ; one finds n ^ while on the 

contrary 0 : 2 (t) t~^. Using the dimensional analysis argument presented above, we immediately conclude that 
q: 2 (t) ~ .j-{ 2 -du- 2 zi^)/(ziy+i) linear quenching of A while for a non-linear ramp a 2 ('r) .^r( 2 -di'- 2 ziz)/{rzi^+i)^ j.j. jg 

noteworthy that whenever vz = 1, the scaling of 0:2 is the same as that of n. This holds true even for the non-linear 

rud 

ramp where both 0:2 and n scale as r . On the other hand, when the environmental XY chain is quenched 

along the gapless Ising critical line h = 1 (varying 7 = t/r across the multi-critical point at h = 1,7 = 0 ), one can 
argue that the multi-critical point plays the role of a dominant critical point with critical exponents that 12 = 1 and 
z = 2m and hence uz 1; substituting these exponents in the KZS relation and the predicted scaling a^., one 
immediately finds n ^ and d ^ . Therefore, the universal scaling formula of d 2 ('/‘) explains all the integrable 

situations discussed in earlier studies. 
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